Entropically-Stabilized Self-compactification in Model Colloidal Systems 
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We discuss the phenomenon of spontaneous self-compactification in a model colloidal system, 
proposed in a recent work on DNA-mediated self-assembly. We focus on the effect of thermal 
fluctuations on the stability of membrane-like self-assembled phase with in-plane square order. 
Surprisingly, the fluctuations are shown to enhance the stability of this quasi-2D phase with respect 
to transition to alternative 3D structures. 
PACS numbers: 64. 70. Kb, 64.70. Nd 



I. INTRODUCTION 
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Effects of physical dimensionality on crystallization are among the most important problems in Condensed Matter 
Physics. As was established by Landau and Pierls pj, the long-range crystalline order in 2D is universally destroyed 
f— — , ' by thermal fluctuations. The problem has been revisited in 1970s by Kosterlitz and Thouless Q-Q, wn0 have shown 
that crystals do exist in 2D, within a new topological definition. Nevertheless, their melting temperature is believed 
to be always lower than in 3D systems (assuming the same interparticle potential). One of the manifestations of this 
^ . effect is the phenomenon of surface melting: a microscopic liquid layer normally appear at the interface of a crystalline 
solid well below its bulk melting temperature. 

In this paper, we describe a remarkable example in which the thermal fluctuation stabilize a 2D crystalline solid, 
embedded in 3D physical space, with respect to transition to an alternative 3D structure. The model system discussed 
*^ | below has been introduced in the recent work by one of us |3J , in order to describe DNA-assisted self-assembly of 
colloids. Its essential ingredients are cohesive interparticle interactions and medium-range soft core repulsion. The 
binary system of same-size spheres (A and B) discussed in Ref. |3j, combines the repulsive potential U(r) acting 
between same-type particles, with A-B attraction. As was shown at that work, both interactions may be induced by 
properly designed DNA. It was found that this colloid-DNA mixture may exhibit an unusually diverse phase diagram 
as a function of two control parameter: the relative strength of attraction and repulsion, and aspect ratio £/a (£ is 
^ ■ the range of repulsive potential U(r), and a is the particle diameter). 

Among various self-assembled phases expected for that system, it was especially striking to find quasi-2D membrane 
with the in-plane square order (SQ). In other words, according to our calculations, this 3D system may prefer to self- 
assemble into a lower-dimensional structure. We will refer to this phenomenon as spontaneous self-compactification. 
Of course, there are other known examples of self-compactified structures in condensed matter, such as lipid mem- 
' branes Q . However, our case is quite unique because it is based on isotropic pair potentials (in contrast to anisotropic 
interactions between lipids, or covalent bonding of carbon atoms in graphite). Note also that at the found SQ phase 
is not a stacking of weakly coupled layers (lamella-like), but rather an isolated membrane-like structure. 

In our early calculations, we have only accounted for the interplay of repulsive and attractive energies, while the 
thermal fluctuations were totally ignored. Even though this approximation is applicable when the characteristic 
energies considerably exceed ksT, the entropic effects are expected to be significant in any realistic case. Given 
that fluctuations are known to strongly affect 2D crystals, one might wonder whether the phenomenon of self- 
compactification will still be present if the fluctuations are introduced. Below we present the detailed study of this 



(j , question. 



II. THE MODEL AND ITS GENERIC FEATURES 



Before going into a specific example, we describe our model and its generic features. We start with particles packed 
into ideal crystalline lattices, whose fluctuation free energies are to be compared. The interparticle potentials will 
be replaced with linear springs, whose spring constants correspond to the second derivatives of the corresponding 
potentials. The first derivatives of the potentials will give rise to a pre-existing stresses in those springs. We consider 
only very short-range interactions, so that any connections beyond second nearest neighbors will be neglected. 

Let k and k be spring constants for the first and second nearest neighbor bonds, respectively. Repulsion between 
the second nearest neighbors induces tension f in K-springs, which should be balanced by an appropriate compressive 
force — r in K-springs. If the interaction range £ is much shorter than nearest neighbors distance a, one expects k»k 
, and r ~ f ~ This gives rise to a hierarchy of elastic constants in this system: k 3> h S> rja. Below, we will 

use harmonic analysis to diagonalize the phonon Hamiltonian and find the fluctational contributions to Free Energies 
of the competing phases. As an example, we compare 2D square lattice (SQ) (embedded in 3D physical space) to an 
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FIG. 1: SQ (a) and its 3D-dual phase (b). 

alternative three-dimensional phase with a very similar local structure. Figure ^shcws this 3D counterpart of SQ 
(referred below as "dual phase"), which also has four nearest neighbors lying in one plane around each particle. The 
difference from SQ is that now there are eight, rather than four second nearest neighbors for each site. Because of 
this difference in the number of repulsive bonds, SQ phase is generally preferred energetically over its 3D dual, at the 
zero temperature limit. Below, we will see whether the free energy balance between the two phases may be reversed 
by the thermal fluctuations. 

Before we present the exact results for SQ, and its 3D-dual phase, we discuss certain general features of this problem. 
Given the hierarchy of the spring constants, one can distinguish between several kinds of phonon modes. Namely, 
there are hard modes, which involve deformations of strong K-springs, and soft modes which only depend on spring 
constant h. There are 3 degrees of freedom per particle, and 2 of them correspond to hard modes, both for 2D lattice 
and its 3D dual phase. In 3D structure, the third mode is the soft one (i«-mode). However, if there were no stresses r 
and f , within harmonic approximation there would be no restoring force for the out-of-plane fluctuations in the 2D 
structure. This means that in 2D lattice the effective spring constant for the transverse modes should be of the order 
of r/a (since i>K> r/a, we may justifiably call these modes of 2D lattice supersoft). Because of the replacement 
of the soft modes with the supersoft ones, we do expect the entropy of the quasi-2D structure to be higher, which will 
further shift the free energy balance to its favor. In addition, in small-g limit, the out-of-plane modes in 2D structure 
will universally have dispersion ~ q A (associated with bending rigidity). This will give another negative correction to 
free energy compared to dispersion ~ q 2 of regular acoustic phonons. Note that these arguments are quite universal, 
and can be applied beyond the particular case of SQ lattice. 

III. EXAMPLE: SQUARE VS. ITS 3D-DUAL 

We now proceed with the discussion of the specific example, 2D SQ vs. its 3D-dual phase. The latter 3D structure 
may be obtained from a simple cubic lattice, by removing the particles occupying sites . (2m, 2k, 21 + 1) and (2m + 
l,2n + 1,2k). In this geometry, the equilibrium condition requires r = —2^2^ , compared to r = —^/2f for SQ. 
Note that it is tension f , not r which is expected to be nearly identical for the two competing structures, since it is 
given by derivative of the potential, —d r U(r) (taken at the distance of the second nearest neighbor). In contrast, r 
can be considerably varied by relatively small deformation of the strong K-spring. 

An obvious choice of unit cell for 2D lattice is a square containing one atom. For the 3D structure, consider a 
cube containing eight smaller cubes, as shown on Figure db). The particles are located at (2m, 21, 2k) + x, where 
x is (1,0,0), (1,0,0), (0,0,1), (1,0,1), (1,1,1), (1,1,0). The translational symmetry in this structure is generated 
by vectors (2,0,0), (0,2,0) and (1,1,1), i. e. the cube on the picture corresponds to two unit cells. Given the 
translational symmetry, the particle at (1, 1, 1) is equivalent to (0, 0, 0) (we label them as type 'a' particles), (0, 1, 1) 
to (1, 0, 0) (type 'b') and (1, 0, 1) to (0, 1, 0) (type 'c'). We can now obtain Hamiltonian H and calculate the phonon 
contribution to free energy. 

The number of particles per a unit cell is, n = 1 for SQ and n — 3 for its 3D-dual. For the 2D phase, displacement 
of a particle at (m,n) is u(m,l) — J^q " (l) e^ mqi+nq2 \ For the 3D phase, the displacements there are three families 
of modes, corresponding to the three types of non-equivalent particles, 'a', 'b', and 'c': 
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u(2m + 1, 2k + I, I) 
u{2m + 1 + 1,2k + 1,1) 
u(2m + l,2k + 1 + 1,1) 




(1) 



The phonon free energy has general form 



= -k B T\ogj D [ui (q)] exp (—^f) (2) 
Here index i = 1, ...,3n parameterizes all the modes for a given wave vector q, and Hamiltonian H is given by 

H = 7« (q) ^ (q) u j (-q) 

q 

By performing the Gaussian integration over Ui (q), one can transform Eq. into: 



= _^ / — 1_ log det 



2?( 



7(q) 



(3) 



The integration here is performed over a single Brilluen zone, i.e. each component of the wave vector runs from — tt 

tO 7T. 

Note that in our regime ( k 3> K and k 3> t/o), one can neglect the coupling between hard and soft modes, which al- 
lows one represent the determinants entering Eq. |[SJ in factorized form: det [7 (q)] ~ det \j( hard ) (q)] det py( so ^*) (q)] . 
Furthermore, the hard modes have essentially identical spectra in the both structures (corresponding to one- 
dimensional chain of K-springs). Therefore, the difference of the fluctuational free energies is mainly due to the 
soft mode contributions. In the case of SQ phase, there is only one soft mode: 



det 



T2D (q)J = 7 ZZ (q) = — — sin J sm 2 J 



For the 3D-dual phase, 3 out of 9 modes are soft, namely (u ax , Ub y ,u c 



4 S D ft) (q)=2«(l + 3a)a« + 



«(l + a) 



x (91,92) X (91,93) 
X (91,92) x(92,93) 
X (91,93) X(92,93) 



(4) 



(5) 



here a = f / (v2fca), x (9, 9') = [1 + cos (9 + 9') — cos (9) — cos (9')], an d 93 = 2g 3 — q 2 — qi. Thus, 
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8(1 + 3a) 



(6) 



Here w x = x (91, 92) + X (91, 93) + X (93, 92), and tu 2 = 2% (91, 92) X 9a) X (93, 92)- Now we can calculate the 
fluctuational contribution into the free energy difference between the two structures: 



AF (/,) = F, 
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(7) 



Consistent with the above general arguments, we obtain F 2D —F^ D < 0, i.e. the entropic effects enhance the stability 
of the self-compactified phase. One can identify the physical origin of each of the three contributions at the above 
result. The first logarithmic term is due to an entropic gain of a single fluctuation particle with all of its neighbors 
fixed. One can see that these fluctuations are enhanced in 2D case (for an arbitrary value of parameter a), mainly 
due to the smaller number of second nearest neighbors there. The second term, —NksT\og2 represent an additional 
gain due to q 4 dispersion of the soft acoustic phonon in SQ phase. In contrast, the soft modes in 3D phase correspond 
to optical phonons, whose dispersion gives rise to the negligible third term in Eq. 0. 
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FIG. 2: 2D honeycomb structure(a) and two of its 3D-dual phases (b), taken from Ref. 



IV. DISCUSSION 



Since the effects responsible for the dominant contributions to our result, Eq. J7J), are not specific for the above 
example, we expect our conclusions to be rather universal, and applicable to other self-compactified phases. In fact, 
the very same binary system introduced in Ref. Q, provides other examples of such phases, e. g. honeycomb lattice. 
This structure has several 3D counterparts [f| , some of which are shown on Fig. © ■ We expect the thermal fluctuations 
to enhance the stability of honeycomb with respect to transition to alternative 3D structures, as in the case of SQ. 

Another aspect of our observation is that the region of the phase diagram corresponding to quasi-2D structures 
(like SQ or honeycomb), should expand when thermal fluctuations are introduced. However, in this study we have 
not discuss the transition between solid and liquid phases. In fact, the strong fluctuations in 2D structures indicate 
that they should melt earlier than 3D solid phases, which is consistent with the classical picture. In other words, the 
thermal fluctuations enhance the stability of quasi-2D latices with respect to transition to 3D crystalline structures, 
but not to the disordered liquid phase. 

Note also that the self-compactified structure need not to be an ideal crystal, even topologically. Introduction of 
defects such as disclination or dislocation, normally associated with 2D melting, is not expected affect our arguments, 
as long as the systems remains effectively two-dimensional. That is because the dominant contribution to the free 
energy is associated with length scales of the order of the interparticle distance. Hence, the macroscopic properties 
and conformation of these lattices are not relevant for the discussed phenomenon. As long at the large-scale behavior 
is concern, one might expect their properties to be similar to those of tethered (solid- like) membranes [fj . It should be 
noted that the non-zero bending rigidity is known to stabilize such a membrane with respect to crumpling transition. 

One may give a simple qualitative interpretation to the obtained results. The self-compactification in zero- 
fluctuation limit corresponds to the regime of the interparticle repulsion which is too weak to destabilize the 2D 
structure, yet strong enough to cause an overall repulsion between two such layers. Introduction of thermal fluctua- 
tions results in an additional effective repulsion between these 2D layers, an effect very similar to Hclfrich interaction 
between conventional membranes , • 
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